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Abstract 

This  document  is  the  final  report  for  AFOSR  Grant  F4962D-99-1-0070,  lattice  Hydrodynamica.**  Under 
the  terms  of  this  grant,  the  Center  far  Computational  Science  at  Boston  University  provided  theoretical 
and  computational  support  to  the  Lattice^Gas  Theory  and  Computation  group  at  the  Space  Vehicles 
Directorate  of  Air  Force  Research  Laboratory  (AFOSR  task  2304CP).  The  principal  research  topics  we 
the  development  of  quantum  lattice-gas  models  for  impienientation  on  quantum  computers,  and  classical 
lattice  Boltzmann  models  for  complex  fluids.  Seven  publications  have  resulted  from  this  effort,  as  well 
as  the  editorship  of  the  Proceedings  of  the  Workshop  on  Quantum  Computation  for  Physical  Modeling 
(QCPM  2000),  which  was  held  in  Falmouth,  Massachusetts  in  October,  2000-  These  proceedings,  edited 
by  Bruce  M.  Boghosian  and  Jeffrey  Yepez,  will  appear  in  the  journal  Computer  Physics  Communications. 
The  principal  results  of  this  work  arc  described  in  this  report. 
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Figure  1:  Two- fermion  probability  density 


1  Principal  Areas  of  Study 

1.1  Quantum  Lattice  Models 

Classical  lattice-gas  automata  (LGA)  are  discrete  dynamical  systems  in  which  a  finite  number  of  state  bits 
are  maintained  at  each  site  of  a  grid,  and  evolve  synchronously  according  to  logical  rules  involving  the  state 
of  a  given  site  and  its  neighbors.  If  the  LGA  has  a  total  of  N  bits  of  state,  it  can  be  in  any  one  of  states. 
The  dynamics  is  thus  a  deterministic  or  stochastic  Markov  process  on  these  states.  Since  the  mid  1980’9,  it 
has  been  known  how  to  construct  LGA  whose  coarse-grained  behavior  is  that  of  a  Navier- Stokes  fluid. 

Quantum  lattice- gas  automata  (QLGA)  are  a  generalisation  of  LGA.  Whereas  an  LGA  with  N  bits 
can  be  in  any  one  of  2^  distinct  states,  a  QLGA  can  be  in  a  complex  superposition  of  all  of  those  2^ 
states  -  although  conservation  of  mass  (and  momentum  if  appropriate)  restricts  the  ewlution  to  particular 
sectors  of  the  full  Hilbert  space.  The  moduli  squared  of  these  complex  amplitudes  may  be  thought  of  as  the 
probability  that  the  system  is  in  the  corresponding  state.  As  these  amplitudes  evolve  according  to  unitary 
transformations,  they  can  inttrfere  with  one  another,  giving  rise  to  behavior  that  is  not  possible  with  classical 
lattice-gas  automata. 

The  AFRL  group  led  by  Jeffrey  Yepez  is  currently  engaged  in  a  program  of  realizing  QLGA  using  nuclear 
magnetic  resonance  (NMR)  apparatus.  The  NMR  apparatus  used  is  very  nearly  “off  the  shelf;  it  required 
very  little  modification  from  those  used  for  medical  imagery.  The  idea  is  that,  by  placing  certain  organic 
molecules  in  a  strong  magnetic  field,  one  can  control  the  spin  States  of  various  nuclear  spins  using  radio 
frequency.  The  nuclear  spins  are  coupled  by  the  electronic  structure  of  the  molecule.  Because  the  NMR 
sample  contains  order  10^®  molecules,  the  measured  signal  effectively  averages  over  these,  yielding  a  high 
signal*to-noise  ratio.  (It  is  to  be  emphasized,  however,  that  only  spins  within  the  same  molecule  enter  into 
a  single  computation.)  The  apparatus  is  able  to  localize  subsets  of  the  molecules  according  to  their  spatial 
position  by  using  field  gradients. 

During  the  granting  period  covered  by  this  report,  Bruce  Boghosian  and  Jeffrey  Yepez  have  coauthored 
a  paper  entitled  “An  Efficient  and  Accurate  Quantum  Lattice-Gas  Model  for  the  Many-Body  Schrddinger 
Wave  Equation,”  which  will  appear  in  the  journal  Computer  Physics  Communications  in  2001  [1].  This  paper 
contains  the  principal  results  obtained  on  the  QLGA  project  durir^g  this  reporting  period;  we  described  and 
simulated  a  quantum  lattice-gas  model  for  a  system  of  particles  governed  by  the  nonrelativistic  Schrodinger 
equation.  The  paper  also  contains  the  first  two-fermion  simulations  of  Schrodinger’s  equation  using  the 
quantum  lattice-gas  algorithm.  An  illustration  of  the  two-partide  probability  density  versus  position  is 
presented  in  Fig.  1. 

In  this  work,  we  also  measured  the  accuracy  of  the  method,  as  a  function  of  spatial  resolution  [l].  We 
found  that  the  method  was  fourth-order  accurate  in  space.  More  significantly,  however,  we  found  that 
the  emergent  Schrodinger  evolution  is  not  exactly  unitary,  as  had  been  previously  reported  [2],  even  for 
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the  simplest  case  of  one  panicle  in  one  dimension.  In  that  case,  the  dynamical  variables  are  and 

and  the  quantity  that  satisfies  the  SchrSdinger  equation  is 

Because  the  method  follows  umtar>-  evolution,  it  is  true  that 


for  t,  but  it  does  not  follow  from  this  that 


t 

In  fact,  the  latter  is  true  only  in  the  continuum  limit.  In  other  words,  while  the  QLGA  is  unitary,  the 
emergent  Schrodinger  equation  is  not. 

These  observations  have  shaped  the  direction  of  our  future  research  on  QLGA.  We  are  currently  adopting 
a  new  model  for  the  QLGA  that  will  yield  perfectly  unitary  Schrodinger  evolution,  by  casting  an  operator- 
splitting  approach  into  the  form  of  a  QLGA. 

Finally,  in  October  of  2000,  Bruce  Boghosian  participated  in  the  organization  of  the  meeting  QCPM  2000 
(Quantum  Computing  for  Physical  Modelling).  Details  about  this  meeting  can  be  found  on  the  meeting  web 
page: 

http :  //xyz .  plh .  af.mil/QCPM2000/QCPK2000_dlab.html 

Jeffrey  Yepe'^  and  Bruce  Boghosian  arc  currently  coediting  the  proceedings  for  this  meeting,  which  will 
appear  in  the  journal  Computer  Physics  Communications.  A  partial  table  of  contents  for  these  proceedings 
appears  as  an  appendix  to  this  report. 

1.2  Entropic  Lattice  Boltzmann  Models 

Over  the  last  ten  years,  lattice  BGK  methods  (a  particular  kind  of  lattice  Boltzmann  model)  has  emerged 
as  one  of  the  more  important  and  successful  methodologies  of  computational  fluid  dynamics.  As  has  been 
noted  by  several  authors,  however,  lattice  BGK  methods  lack  an  if-theorem,  and  are  plagued  by  a  variety 
of  Dl-understood  numerical  instabilities.  In  recent  work,  we  have  demonstrated  a  generalization  of  the  BGK 
techni'^ue  that  makes  it  possible  to  reintroduce  an  //-theorem  in  many  cases.  During  the  period  of  this 
granv,  v.-^e  have  completed  and  published  a  general  program  for  the  construction  of  “cntropically  stabilized’’ 
ht'Le  Boltzmann  models,  and  illustrated  its  application  to  several  example  problems  [3]. 

The  crux  of  the  idea  is  to  encourage  the  model  builder  to  specify  an  appropriate  H  function  (Lyapunc^’ 
function)  for  the  model,  rather  than  tr>'  to  blindly  dictate  an  equilibrium  state,  as  is  done  for  BGK  models. 
Of  course,  specifying  a  form  for  M  determines  the  equilibrium  distribution,  but  it  also  governs  the  approach 
to  this  equilibrium.  It  can  therefore  be  used  to  control  the  stability  properties  of  the  model.  It  should  be 
emphasized  that  the  presence  of  a  Lyapunov  function  guarantees  the  nonlinear  stability  of  the  model,  which 
is  a  much  stronger  condition  than  linear  stability. 

In  earlier  work,  under  AFOSR  grant  number  F49620-95-1-0285,  we  showed  that  collision  operators  that 
are  constructed  to  increase  H  defined  in  this  way  are  similar  in  form  to  the  lattice  BGK  collision  operators, 
except  that  their  relaxation  parameter  maj’  depend  on  the  current  state.  As  a  result,  the  transport  coefficients 
may  have  a  certain  minimum  value  in  models  of  this  type  but  these  minimum  values  are  often  actually  zero. 
The  results  from  that  study  are  illustrated  in  Fig.  2,  which  shoe's  the  decay  of  the  amplitude  of  an  initially 
sinusoidal  density  profile  in  a  simple  entropic  lattice  Boltzmann  model  of  diffusion.  T^ble  1  then  shows  the 
theoretical  value  of  the  diffusion  versus  the  value  measured  from  numerical  experiments.  We  see  that  the 
agreement  is  good  for  values  of  the  diffusivity  (in  lattice  units)  greater  than  one  part  in  10^. 

In  the  present  work,  we  have  demonstrated  that  this  makes  possible  the  construction  of  fully  explicit, 
perfectly  conservative,  absolutely  stable  algorithms  with  arbitrarily  small  transport  coefficients.  In  fact,  we 
have  showed  that  the  ultimate  limitation  to  the  application  of  these  algorithms  for  very  small  transport  co¬ 
efficients  come  from  considerations  of  accuracy,  rather  than  stability.  Wg  have  created  cntropically  stabilized 
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Figure  2:  Decay  of  ainusoida]  density  profile  for  k  =  0.99999 
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Figure  3:  Decay  of  ainufioidal  density  proflJe  for  k  =  0-99999 
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T^ble  1:  Theoretical  and  meaanred  values  of  the  diffiisivity  for  various  values  of  k. 
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lattice  Boltzmann  models  of  two  dimensional  Navier-Stokes  flow  that  arc  stable  for  values  of  the  viscosity 
in  lattice  units  that  are  more  than  two  orders  of  magnitude  smaller  than  what  is  possible  in  lattice  BGK 
simulations.  Since  computational  complexity  scales  as  Reynolds  number  cubed,  this  implies  s  six-order  of 
magnitude  decrease  in  the  number  of  site  updates  necessary.  The  entropically  stabilized  algorithm  requires 
more  than  two  orders  of  magnitude  more  floating-point  operations  per  site  update.  Kevertheless,  the  net 
result  would  be  a  three-order  of  magnitude  improvement  in  computational  complexity  for  a  given  calculation 
-  if  stability  were  the  only  limiting  factor.  In  practice,  other  accuracy  considerations  may  enter;  on  the  other 
hand,  the  entropically  stabilized  algorithm  may  well  be  valid  ’.-eyond  the  point  where  the  smallest  eddies 
are  on  the  order  of  the  grid  size.  In  some  sense,  cntropic  stabilization  may  constitute  a  new  kind  of  subgrid 
turbulence  model,  and  this  is  driving  the  direction  of  pres'=^r,ii  research  on  this  topic.  Our  publication  gives 
more  details  on  this  issue  [3]. 

In  order  to  conduct  these  studies,  we  had  to  develop  algorithms  for  cfBcient  site  updates  using  the 
entropically  stabilized  algorithm.  There  are  two  principal  algorithmic  problems  that  enter  for  models  of 
fluids  in  two  and  three  dimensions.  First,  it  is  first  necessary  to  find  the  equilibrium  distribution  for  given 
values  of  the  conserved  quantities.  Second,  it  is  necessar)'  to  find  the  minimum  value  of  the  collision  parameter 
Train  foT  which  the  step  will  be  stable.  The  value  of  the  collision  parameter  r  actually  used  is  then  rniin/«i 
where  /c  <  1.  We  have  shown  that  the  transport  cocfflcicnt  goes  to  zero  as  «  I- 

To  solve  the  first  of  these  problems,  we  noted  that  the  equilibrium  distribution  could  be  perturba- 
tiveb'  expanded  in  Mach  number  [3],  Subsequent  to  the  publication  of  our  paper,  we  found  that  if  this 
Mach-number-expanded  equilibrium  is  used  as  the  initial  guess  in  a  Newton-Raphson  solution  for  the  exact 
equilibrium,  the  solution  will  converge  in  one  pass  of  the  Newton-Raphson  algorithm  for  Mach  numbers  less 
than  0.3. 

At  the  time  of  publication  of  the  paper  [3],  we  had  no  simQar  Mach  number  expansion  of  Train-  In  subse¬ 
quent  work,  we  found  that  the  difficulty  stemmed  from  the  nonanalyticity  of  Tmin  at  the  equilibrium  point. 
This  is  illustrated  for  the  diffusion  model  in  Fig.  4.  We  were  then  able  to  construct  rational  approximations  - 
Pad6  approximations  -  to  Tmin  that  could  likewise  be  used  as  the  initial  guess  in  a  Newton-Raphson  iteration 
for  the  exact  solution.  As  with  the  equilibrium  point,  we  have  found  that  such  Newton-Raphson  solutions 
converge  in  one  step,  thereby  making  possible  efficient  entropic  lattice-Boltzmann  simulations.  These  latter 
algorithmic  advances  will  be  reported  on  in  a  future  publication  [4]. 

2  Conclusions 

We  have  described  the  program  of  activity  of  the  lattice-gas  research  collaboration  between  the  Boston 
University  Center  for  Computational  Science  and  the  Space  W\icles  Directorate  of  AFRL  (under  AFOSR 
task  2304 CP).  This  collaboration  has  centered  on  three  major  areas  of  study:  (i)  lattice-fluid  models  for 
droplet  formation,  (ii)  entropic  lattice  Boltzmann  models,  auo  (iii)  quantum  lattice-gas  automata.  We  have 
provided  a  detailed  account  of  the  principal  new  results  in  all  three  of  these  areas.  Nine  publications  have 
resulted  from  this  effort,  as  well  as  the  sponsorship  of  the  Seventh  International  Conference  on  Discrete 
Models  for  Fluid  Mechanics,  and  the  preparation  of  the  proceedings  of  that  meeting  as  a  special  issue  of  the 
International  Journal  of  Modern  Physics  C. 
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Figure  4:  Minimum  tau  as  functicai  in  allowed  polytope  of  states.  Note  nonanalyticity  at  equilibrium  point. 
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A  Publications 

The  foUcrwing  seven  publications  coauthored  by  Bruce  M.  Boghosian  were  related  to  the  topic  of  this  grant 
proposal,  and  were  published  during  the  time  period  of  the  grant: 

1.  B.M.  Boghosian,  J.  Yepez,  P-V.  Coveney,  AJ.  Wagner,  “En tropic  Lattice  Boltzmann  Methods Proc. 
Roy.  Soc.  Lon.  A  457  (2001)  717-766-  coii<i-inat/0005260 

2.  M-  Nekovee,  P.V.  Coveney,  H.  Chen,  B.M.  Boghoeian,  “A  Lattice-Boltzmaiin  Model  for  Interacting 
Amphiphilic  Fluids,”  Phys.  Rev.  E  62  (2000)  8282-6294.  cond-Bat/0006319 

3.  Francis  J.  Alexander,  Bruce  M.  Boghoeian,  Richard  C.  Brower,  S.  Roy  Kimura,  “Fourier  Acceleration 
of  Langevin  Molecular  Dynamics,”  to  appear  Pkys.  Rev.  E  (2001).  cond-mat/0001418 

4-  B-M,  Boghosian,  C.  Chow,  T.  Hwa,  *Hydrodynamics  of  the  KuramotoSivashinsky  Equation  in  Two 
Dimensions,”  Phys.  Rev.  Lett.  83  (1999)  5262-5265-  coad-aat/99 11069 

5.  B.M.  Boghosian,  P.V.  Coveney,  “A  Particulate  Basis  for  an  Immiscible  Lattice-Gas  Model,”  Comp. 
Phys.  Comm.  129  (2000)  46-55.  cond-ciat/9911340 

6.  H.  Chen,  B-M.  Boghosian,  P.V.  Coveney,  M.  Nekovee,  “A  Lattice  Boltzmann  Model  of  Ibmary  Am¬ 
phiphilic  Fluids,”  Proc.  Roy.  Soc.  London  A  450  (2000)  2043-2057.  cond-mat/99 10369 

7.  B.M.  Boghosian,  P.V.  Coveney,  P.J.  Love,  “A  Three-Dimensional  Lattice-Gas  Model  for  Amphiphilic 
Fluid  Dynamics,”  Proc.  Roy.  Soc.  Lon.  A  450  (2000)  1431-1454.  cond-mat/99 07 298 
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B  Invited  Talks  and  Presentations 

The  following  invited  tallcs  and  presentations  by  Bruce  M-  Boghosian  were  related  to  the  topic  of  this  grant 
proposal,  and  were  presented  during  the  time  period  of  the  grant: 

1.  Speaker,  ‘'Symposium  Celebrating  the  Tenth  Anniversary  of  the  Center  for  Computational  Science,” 
Boston  University  (17  November,  2000). 

2.  Speaker,  Greater  Boston  Statistical  Mechanics  Meeting,  Brandeis  Univ'ersity  (14  October  2000). 

3.  Member  of  Scientific  Program  Committee,  "^inth  International  Conference  on  Discrete  Models  for 
Fluid  Dynamics”  Santa  Fe,  New  Mexico  (21-25  August,  2000). 

4.  Speaker,  Gordon  Conference  on  Computational  Physics  Education,  Plymouth,  New  Hampshire  (11-16 
June,  2000). 

5.  Speaker,  Materials  and  Manufacturing  Directorate,  Air  Fbrce  Research  Laboratory,  Wright*Patterson 
Air  Force  Base,  Ohio  (30  May,  2000). 

6.  Minisymposium  Speaker,  “Role  of  Curvature  and  Hydrodynamics  in  Soft  and  Biological  Matter,”  Third 
SIAM  Conference  on  Mathematical  Aspects  of  Materials  Science,  Philadelphia,  Pennsylvania  (22-24 
May,  2000). 

7.  Colloquium  speaker,  Department  of  Physics,  Clark  University,  Worcester,  Massachusetts  (27  April, 
2000). 

S.  Colloquium  speaker,  Department  of  Physics,  Boston  University,  Boston,  Massachusetts  (4  April,  2000). 

9.  Speaker,  “Recent  Developments  in  Computer  Simulation  Studies  in  Condensed  Matter  Physics,’’  Cen¬ 
ter  for  Simulational  Physics,  Department  of  Physics  and  Astronomy,  University  of  Georgia,  Athens, 
Georgia  (21-25  February,  2000). 

10.  Colloquium  speaker,  Department  of  Mathematics,  Tufts  University,  Medford,  Massachusetts  (17  Febru¬ 
ary,  2000). 

11.  Colloquium  speaker,  Department  of  Mathematics,  Worcester  Polytechnic  Institute,  Worcester,  Mas¬ 
sachusetts  (3  February,  2000). 

1/  Seminar  speaker,  Department  of  Applied  Mathematics,  University  of  North  Carolina,  Chapel  Hill, 
North  Carolina  (6  November,  1999). 

13.  Seminar  speaker,  Department  of  Chemistry,  University  of  London,  Queen  Mary  and  Westfield  College, 
London,  United  Kingdom  (9  July,  1999). 

14-  Seminar  speaker,  Center  for  Nonlinear  Studies,  Los  Alamos  National  Laboratory,  Los  Alamos,  New 
Mexico,  10  June  1999. 

15.  Organizing  Committee  member,  NIST  CTCMS  Workshop  on  Hybrid  Methods  in  Multiscale  Modeling 
of  Materials,  12-14  May  1999. 

16.  Speaker,  Molecular  Simulation  ’99,  Advances  in  Mesoscale  Simulation  Methodology,  Web-based  real¬ 
time  conference  at  URL  http://mol«ia.vei.co.\ik/  (27  April  1999). 

17.  Seminar  speaker,  Department  of  Mechanical  Engineering,  Massachusetts  Institute  of  Technology  (6 
AprU  1999). 
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C  Honors  and  Awards 

Bruce  Boghosian  was  elected  a  Fellow  of  the  American  Physical  Society  in  2000.  He  was  also  invited  to  join 
the  editorial  advisory  board  of  the  journal  Physiea  A.  He  remains  on  the  editorial  boards  of  the  International 
JovmaJ  of  Modem  Physics  C,  and  of  Computers  in  Science  and  Engineering. 


D  Proceedings  of  the  Workshop  on  Quantum  Computation  for 
Physical  Modeling  (QCPM  2000) 

This  is  a  list  of  the  publications  that  will  appear  in  the  Proceedings  of  the  Workshop  on  Quantum  Com¬ 
putation  for  Physical  Modeling  (QCPM  2000),  held  in  Falmouth,  Massachusetts  in  October,  2000.  These 
proceedings,  edited  by  Bruce  M.  Boghosian  and  Jeffrey  Yepez,  will  ^peax  in  the  journal  Computer  Physics 
Communications. 


CONTENTS 


Lattice  Bolumann  Schemes  for  Quantum  Applications 
5.  Sued 


A  Quantum  Full  Adder  for  a  Scalable  Nuclear  Spin  Quantum  Computer 
G.P,  Berman,  G.D,  Doolent  G.V.  Ldpez,  V.L  Tsifrinovich 


Quantum  Computing  with  Spin  Qubits  in  Semiconductor  Structures 
y.  Pnvman,  D.  Mozyrsky,  LD.  Vagner 


Simulating  Fermions  on  a  Quantum  Compute 
G  Jrtiz,  J,E.  Guhernatis,  E.  KniU,  R,  Laflamme 
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